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ABSTRACT 


/ 

Internal  stress  in  and  around  an  oblate  spheroidal  inclusion  is  evaluated 
by  using  Eshelby  theory.  We  consider  three  different  sources  of  internal 
stress  as  a function  of  the  aspect  ratio  and  elastic  moduli  ratio;  these  are 
(1)  misfit  effect  (2)  inhomogeneity  effect  and  (3)  plastic  deformation  effect. 
The  misfit  effect  arises  from  the  difference  of  thermal  expansion  coefficients 
of  the  matrix  and  inclusion.  Normal  stresses  within  the  inclusion  (or  the 
principal  stresses),  maximum  shear  stress  inside  the  inclusion,  total  strain 
energy  and  the  normal  stress  at  the  matrix-inclusion  boundary  are  determined. 
The  inhomogeneity  effect  arises  from  the  difference  of  the  elastic  moduli  of 
the  matrix  and  inclusion,  which  perturbs  otherwise  uniform  applied  stress.  The 
direction  of  applied  stres^S  is  parallel  to  the  axis  of  revolution  of  the  in- 
elusion.  Normal  stresses  within  the  inclusion  and  at  the  matrix-inclusion 
boundary  are  obtained.  Plastic  deformation  of  the  matrix  in  the  presence  of  a 
plastically  non-deformable  inclusion  produces  internal  stresses,  as  shown  by 
Ashby  and  by  Tanaka  and  Mori.  The  approach  developed  by  Tanaka  and  Mori  is 
extended  here  to  evaluate  the  internal  stresses  for  the  inclusion  geometries 
appropriate  to  the  analysis  of  technically  important  cases. 


The  results  of  the  present  calculations  in  limiting  cases  agree  with 
previously  published  calculations.  Implications  of  the  results  are  dicussed 
in  connection  with  several  experimental  studies. 


1 . INTRODUCTION 


Inclusions  in  a material  are  known  to  exert  significant  influence  on 
its  mechanical  properties,  especially  ductility  and  fracture  strength.  Typical 
inclusions  have  different  physical  properties,  such  as  elastic  constants  and 
thermal  expansion  coefficient,  from  those  of  the  matrix,  and  internal  stress 
develops  around  the  inclusions  upon  stressing  or  temperature  changes.  Important 
in  the  study  of  the  mechanical  properties  of  inclusion  bearing  materials  are 
the  following  three  internal  stress  problems  caused  by  the  presence  of  an 
inclusion:  (1)  Misfit  Effect:  How  high  internal  stresses  can  be  produced  in 

and  around  the  inclusion  which  has  a misfit  strain?  (2)  Inhomoqeneity  Effect: 

How  does  the  inclusion  perturb  an  external  stress  applied  to  the  material? 

(3)  Plastic  Deformation  Effect:  If  the  inclusion  is  assumed  to  be  deformed 

only  elastically,  how  much  internal  stresses  can  be  developed  in  and  around  the 
inclusion  when  the  matrix  is  deformed  plastically?  The  most  common  example  of 
the  misfit  effect  is  one  due  to  a difference  in  the  thermal  expansion  coefficients 
of  matrix  and  inclusion.  The  case  of  a spherical  inclusion  was  thoroughly  studied 
by  Mott  and  Nabarro  [1].  The  inhomogeneity  effect  has  been  investigated  in 
order  to  obtain  a stress  concentration  factor.  A penny-shaped  cavity  and  an 
elliptical  hole  were  treated  initially  and  have  been  most  widely  considered  [2,3]. 
Sadowsky  and  Sternberg  [4]  extended  the  treatment  to  a prolate  spheroidal  cavity, 
and  Nauber  [5]  obtained  the  solution  for  a spheroidal  cavity.  Edwards  [6]  solved 
the  misfit  and  inhomogeneity  problems  for  a prolate  spheroidal  inclusion. 

Recently,  Mirandy  and  Paul  [7]  evaluated  the  inhomogeneity  effect  for  an  ellip- 
soidal cavity  by  using  Eshelby  theory  [8].  The  Eshelby  theory  of  the  elastic 
strain  field  in  a constrained  transformation  also  was  the  basis  of  the  Tanaka  and 
Mori  analysis  of  the  plastic  deformation  effect  [9].  Their  method  was  applied 
to  the  analysis  of  inclusion  fracture  and  decohesion  in  the  limiting  case  of  the 
inclusion's  shape  [10,11].  Additional  studies  of  this  effect  as  it  pertains  to 


work  hardening  of  multiphase  materials  have  also  ber"  published  [12,13]. 

In  surveying  previous  works  on  the  internal  stress  problems  due  to  an 
inclusion,  it  has  become  clear  that  no  detailed  study  has  been  made  for  an 
oblate  spheroidal  inclusion.  Yet,  this  is  central  in  importance  in  the 
applications  involving  the  ductility  and  fracture  strength  of  a material. 
Therefore,  we  wish  to  evaluate  for  an  oblate  spheroidal  inclusion  with 
respect  to  the  misfit,  inhomogeneity,  and  plastic  deformation  effects. 

Eshelby  theory  [8]  is  used  to  obtain  the  internal  stress  inside  the  inclusion 
and  at  the  matrix-inclusion  boundary.  Effects  of  the  aspect  ratio  of  the 
inclusion  and  differences  in  the  elastic  moduli  of  the  matrix  and  inclusion 
are  examined.  When  applicable,  the  results  obtained  are  compared  with  the 
previously  published  one.  Also  described  are  applications  to  the  analysis 
of  thermally  induced  dislocation-punching  from  a spherical  inclusion,  and  to 
the  internal  stress  evaluation  of  a structural  steel  containing  MnS  inclusions. 

2.  INCLUSION 

Let  X.  (i  = 1,  2,  and  3)  be  Cartesian  coordinate  axes.  Consider  an 
oblate  spheroidal  inclusion  at  the  origin  (see  Fig.  1).  Its  volume  is  Vj, 
and  its  boundary  is  given  as, 

(X2  + X2)/a2  + X2/c2  = 1 . (1) 

The  shape  of  an  inclusion  can  be  described  by  the  aspect  ratio  k(=c/a).  The 
elastic  properties  of  the  inclusion  and  matrix  are  assumed  to  be  isotropic. 
Elastic  stiffness  of  the  matrix  is  given  by  Cij-k£  = AS^S^  + p(6ik6j£  + <5^*$^) 
and  that  of  the  inclusion  by  C^^  = X*<5ij-6k£  + U*(<5ik6j£  + k5iJZ.^  where  * and 
A*  are  the  Lam6  constants,  p and  p*  are  the  shear  moduli,  and  6..  is  the 

' J 

Kronecker  delta.  In  subsequent  sections,  the  Young's  modulus  (E  or  E*)  and 
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Poisson's  ratio  ( v or  v*)  will  also  be  used  because  of  convenience. 


3.  INTERNAL  STRESS;  ESHELBY  THEORY 

Eshelby  [8]  developed  a method  to  analyze  the  elastic  field  in  and 
around  an  ellipsoidal  inclusion,  within  which  exists  a uniform  stress-free 

T* 

strain,  e. ..  In  this  study  we  use  this  theory  to  evaluate  the  internal 
* j 

stresses  of  the  inclusion  and  matrix. 


3-1.  Internal  Stress  of  Inclusion, 

Eshelby  showed  that  the  internal  stress  of  an  inclusion  can  be  written 
as. 


I _ r / c T,  r * , c T*, 

CTi j " Sjkrek£  ' ®k2;  Lijkrek£  " ek£  ' 


(2) 


where  is  total  (or  constraint)  strain  and  ej^  is  stress-free  strain  of  the 

equivalent  inclusion  [8].  In  Eq.  (2)  and  in  what  follows,  the  summation  of 

★ 

repeated  indices  is  implied.  In  the  limiting  case  of  C1jkJt  = CijkH’  we  have 


k£ 


#T* 

®k«. 


T* 

, hence  we  obtain  immediately  for  a given  with  the  use  of 


Eshelby  tensor,  as 


'U 


Sk£mn  emn 


(3) 


In  general,  however,  t cijk«,»  and  n is  necessary  to  obtain  eij  for  a 

T* 
ij 


T*  T* 

given  e. . using  Eq.  (2).  When  the  off-diagonal  components  of  e,.  are  zero 
1 J * u 


and  e 


11 


T*  T 

*22’  eij  becomes  simply 


•1, 


eT„  . « 

22 


'33 


= z, 


(4) 


and  all  the  off-diagonal  components  vanish  Here,  the  equality  ej-j  = 
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results  from  the  symmetry  of  the  inclusion  and  e. Therefore,  remaining 

1 *3 

tasks  in  the  calculation  of  aj,  are  to  obtain  (x,z)  as  a function  of  the 

* J 

T* 

ratios  of  elastic  moduli;  i.e.,  p*/p  and  v*/v  and  k as  well  as  e- .. 

■ J 
M 

3-2.  Internal  Stress  of  the  Matrix,  g}. 

* J 

As  the  distance  from  the  matrix-inclusion  boundary  increases,  the 


internal  stress  of  the  matrix  decreases.  Thus,  the  stress  state  of  the  matrix 
at  the  boundary  is  of  interest.  Because  of  the  symmetry  about  x^  axis,  let 
us  consider  o?.  on  the  plane,  x.  = 0.  At  point  D on  the  boundary  (see  Fig.  1), 

* J 

the  normal  rT,  is  given  by 


n = (cose,  0,  sine). 


where  0 is  the  angle  between  the  ^ axis  and  the  tangential  line  passing 
through  the  point. 

Eshelby  [8]  deduced  the  total  strain  of  matrix  just  outside  inclusion, 

M c T — 

e. .,  which  remains  to  be  elastic,  as  a function  of  e. • , e. ■ and  n.  It  follows 

1 J ■ J ' J 

that  the  components  of  the  internal  stress  of  the  matrix  at  the  boundary  on 


the  *2  Plane  are  given  by 


= u(2(3e^  + «33)  - (2  + 7 cos2e  - 3 cos4e)x  - (2  - 3 cos2e  f 3 cos4e)z}i 


?22  ■ y{2(3e^i  + e!^)  - (2  + cos2e)x  - (2  - cos2e)z} 


M = ..fJ/.C  i 


p{4(e^  + ) - (4  - 4 cos  0 + 3 cos  e)x  - (4  - 3 cos  0)z) 


M _ M _ n 
°12  ‘ °23  ‘ 0 


M2  ? 

o13  = y{(3  cos  0 - 4) X - (3cos20)z)  sine 


cose  . 
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4.  MISFIT  EFFECT 


Let  us  consider  the  case  in  which  an  inclusion  is  subjected  to  a 
hydro-static  misfit  strain,  eT6...  It  can  be  written  as 

* J 


T*  T T 

ei j = C 5ij  0 0 

f 

0 eT  0 


^0  0 e 


4-1.  Misfit  Induced  Internal  Stress  Within  an  Inclusion, 

T*  T T i J i s 

Substituting  e—  = c 6^  into  Eq.  (2),  due  to  misfit  are  obtained 

as. 


m{(l  - m)h(k)  + 1} J 

{(1  - m)2f(k)  - (1  - m)g(k)  + 1} 


m{(l  - m)K,(k)  + 1 } 


{(1  - m)  t(k)  - (1  - m)g(k)  + 1) 


(7) 


where  m is  b*/u,  and  f(k)  = + s-| 1 22^S3333  " 2S1 1 33  S3311 

g(k)  = Snil  + Sn22  + S3333, 

= S1 1 33  " S3333’ 
t(k)  = 2S33n  - Sni]  - S1]22. 


For  k = 0 or  1 , these  terms  of  the  Eshelby  tensor  are  given  elsewhere 
[13,14].  Integration  of  the  original  expressions  is  necessary  for  other 
values  of  k.  The  values  of  f,  g,  h,  and  l have  been  evaluated  for  k = 0.01 


- 5 - 


to  1 and  are  presented  in  Fig.  2.  In  order  to  simplify  the  calculations, 

it  is  assumed  that  v = v*  = 1/3  holds. 

From  the  above  results,  the  internal  stress  of  the  inclusion,  (c^.)  . , 

v lj'mis’ 

can  be  obtained  as  follows: 


(oj , ) • = (olo ) • = 
v 11  mis  ' 22 'mis 


4^ — v-  • {-4e;T  + (3x  . + z . )} 

[1  - m)  v mis  mis' 


(o33>m1s  ’ 2(l3-*m)  -<-^T  * (xmis  + W> 


(o1J>nris  * 0 "hen  1 * J • 


Note  tnat  no  singularity  exists  at  m = 1,  because  the  (1  - m)  term  factors 

out  upon  substitution  of  x and  z . into  Eq.  (8). 

mis  mis 

Results  of  the  calculations  are  presented  in  Figures  3 and  4.  Stresses 
are  expressed  in  dimensionaless  forms  by  dividing  them  by  E*e'r  or  p*4.  The 
variations  of  (aj^mis  and  (032pmis  a9ainst  k are  shown  in  Fig*  3 for  three 
values  of  m.  One  can  see  that,  as  the  inclusion  approaches  a thin  disc 
shape,  the  component  normal  to  the  plane  of  the  disc,  (033)^5,  approaches 
zero,  and  (oj])m1-s  and  ^°22^mis  aPProacl1  an  asympototic  value  of-1.5  e^E* 
regardless  of  m.  The  magnitude  of  dimensionless  stresses,  (0 J -j  )mi s/E*cT 
and  (o^)mis/E*eT,  becomes  larger  as  m decreases.  The  maximum  shear  stress 
in  the  inclusion  is  given  by 


Tmis  ~ ^a33^mis  " ^°l4mis}/2 


— 4 ■ 

= p*  • e 


mU(k)  - h(k)} 

{(1  - m)  f (k)  - (1  - m)g(k)  + 1} 


IT  I 

Figure  4 show  plots  of  x t /p*e  against  k.  At  k = 1 , x!  vanishes  due  to 
3 mis  mis 

the  equality  (oJ,)m1s  ■ (o‘2)mjs  - (o‘3)mfs.  «hen  k becomes  zero,  tJs 

reaches  regardless  of  m.  Note  that  the  effect  of  m is  smaller  on 

x*.  than  on  the  normal  stress  components, 
mis 

4-2.  Strain  Energy  Due  to  Misfit  Strain 

Once  the  diagonal  components  of  (a|j)mi-s  are  known,  it  is  easy  to 

obtain  the  strain  energy  E$  due  to  the  misfit  effect.  E$  can  be  written  simply 

as,  E /VT  - - 0.5  eT(o.  J).  . When  neither  p*/y  nor  v*/v  is  unity, 

Si  kk  mi s 

Es/Vj  = p*(eT)2  • P/Q 

where 

P = (1  + v)(l  + v*)[3(l  -v)  + 3(1  - m)  (1  - v)  {f (k)  - g ( k ) } - m(l  + v)] 
and 

Q = [(1  - m)  (1  + v)(l  - v)  ( 1 - 2v*)  {f (k)  - g(k) ) - m(l  - m)(l  - 2v)(l  + v*)  f (k) 

+ (1  + v){m(v*  - v)  + (1  - v)  (1  - 2v*)}]. 

At  k = 0,  we  have  f (0)  = 0 and  g(0)  = 1.  Thus,  P and  Q are  simplified  as 

P = 2(1  + v)(l  + v*)(l  - 2v)m 

and 

Q = (1  + v)(l  - v*)(l  - 2v)m, 

which  yield 

Es/Vi = 2u*  ' * (fJ)2-  (10^ 
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This  result  is  apparently  equal  to  that  derived  by  Robinson  [15],  except 
that  the  Poisson's  ratio  in  the  numerator  is  that  of  the  matrix  in  his 
expression.  However,  this  discrepancy  appears  to  be  a typographical  error. 

Regardless  of  the  aspect  ratio,  when  and  p are  equal,  can  be 
given  by  Eq.  (10).  E$  is  independent  of  the  shape  of  the  inclusion,  as 
originally  deduced  by  Eshelby  [8].  In  his  derivation,  however,  the  Poisson's 
ratios  of  the  inclusion  and  matrix  were  taken  to  be  identical.  For  more 
general  cases,  aspect  ratio  dependence  of  Es  is  presented  in  Fig.  5,  taking 
both  v and  v*  to  be  1/3.  Whereas  Eg  is  independent  of  k when  m is  unity,  it 
is  interesting  to  note  that  a spherical  inclusion  has  a minimum  E$  for  m > 1 
or  a maximum  E$  for  m < 1.  In  agreement  with  findings  of  Barnett  et  al . [16], 
E$  is  larger  in  discs  than  in  a sphere  when  m > 1 (or  u*  > u). 

4-3.  Internal  Stress  of  the  Matrix 
4-3-1.  Equator  and  Polar  Points 

The  internal  stress  of  the  matrix  at  the  matrix-inclusion  boundary 

M 

varies  from  point  to  point.  Here,  we  present  numerical  values  of  (033)mis 

at  an  equator  point  A and  a polar  point  B (see  Fig.  1).  By  using  x^. ? and 

z . given  in  Eq.  (7),  (oM.)  , at  A and  B can  be  obtained  from  Eq.  (6). 
mis  l j mis 

The  results  for  (cu-J  • are  shown  in  Fig.  6.  Note  that  at  the  polar  point 

j j ni  i b 

•j  t |y| 

B,  (°33)mis  is  equal  to  (033)^  (cf.  Fig.  2).  As  k becomes  smaller,  (a33)mis 
depends  more  on  m at  the  equator  point  A.  At  the  polar  point  B,  effect  of  m 

M 

on  (o„)  •„  increases  as  k approaches  unity. 

33  mis 

Edwards  [6]  calculated  (0^.)  . at  A and  B for  prolate  spheroids;  i.e., 

ij  mis 

M M 

k >_  1 . He  obtained  (o^-])m^s  at  A and  (o33)m^s  at  B.  These  are  both  the 
normal  component  at  the  boundary,  and  can  be  obtained  from  • easily.  Our 
results  agree  with  his  results  at  k = 1,  except  for  small  discrepancies  due 
to  the  different  values  of  v.  In  our  calculations,  v = v*  = 1/3  is  used 
whereas  Edwards  took  v = v*  = 0.3. 
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4 - 3 - 2 . Maximum  Shear  Stress  of  the  Matrix 

When  an  inclusion  takes  spherical  shape,  (k  = 1),  xmi$  and  z^s  in 
Eq.  (7)  become  identical  because  of  an  equality  h(l)  = £.(!'.  Hence,  (°^j)m-js 
at  the  matrix-inclusion  boundary  on  the  plane  x^  = 0 becomes 


(a^-i)  . = |'P(1  ' cos20)x  . 

v ll'mis  3 ' ' mis 


^°22^mis  3*P  Xmis 


^a33^mis  = ' ^*y(2  " 3 cos20)xmis 


(ID 


, M » _ i M , 

(o12'mis  " (a23)mis 


^°13^mis  = - 4y  sin6C0S9  xmis  • 


The  principal  stresses  of  Eq.  (11)  are  given  by  a rotation  about  the 
x2  axis  by  (-6)  as 


M -10  . ,,  y 
°1  = ~T  y Xmis 


oM  . 4 . x . 
&2  " 3 1 mis 


^ * ! * V1  xmis 


and  the  maximum  shear  stress  of  the  matrix  at  the  boundary  becomes 


max1 


1 , M M 

2 I°1  ” °3 


^Vis1 


(12) 


9 - 


M 

This  corresponds  to  ( 0^3 )m-j s in  Eq.  (11)  with  0 = 45°.  At  k = 1,  in 
Eq.  (7)  becomes 


x - 3m  T 

mis  (2m  + TJ'  * e ’ 


(13) 


hence,  we  obtain 


max' 


6m 

1 + 2m  M 


(14) 


The  factor  6m/(l  + 2m)  increases  from  0 to  2 for  a change  in  m from  0 to  1 . It 
gradually  increases  to  the  asymptotic  value  of  3 as  m diverges. 

This  particular  problem  was  treated  first  by  Mott  and  Nabarro  [1]  via 
the  method  of  "a  centre  of  dilatation"  [17,18].  Suppose  an  elastic  sphere 
of  a radius  (1  + eT)2Q  is  forced  into  a spherical  hole  of  a radius  aQ  in  an 
infinite  matrix.  After  removing  surface  tractions,  which  squeeze  the 
sphere  to  radius  aQ,  its  equilibrium  boundary  relaxes  to  a radius  a = (1  + e)aQ. 
Let  K*  be  the  bulk  modulus  of  the  sphere.  Also  "the  effective  bulk  modulus 
for  expanding  a hole"  is  taken  to  be  4p/3  [19].  From  the  balance  of  a 
radial  stress,  we  have  as  the  first  approximation, 


3K* 

‘ 3K*  + 4p 


(15) 


The  stress  at  outside  the  sphere  is  given  by 


“tj  ■ 2 “ • ‘ • *3<sij/r3  - 3xfVr5) 

where 

r = (X2  + X2  + Xj)172. 


(16) 
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The  maximum  shear  stress  is  obtained  at  (a//2  , 0,  a//2  ),  and  is  given  by 

M _ -9K*  T 

Tmax  " 3K*  + lil  * y * e ‘ (17) 

When  v*  = 1/3,  Eq.  (17)  reduces  to  Eq.  (14). 

Let  us  apply  Eq.  (14)  to  the  analysis  of  dislocation  generation  due  to 
thermal  stresses.  Consider  that  a spherical  inclusion  has  a thermal  expansion 
coefficient,  a*,  different  from  that  of  the  matrix,  a.  The  inclusion  is 
subjected  to  a stress-free  transformation  strain  e*"<5. . with  = (a*  - a)AT 

* J 

after  cooling  from  a high  temperature.  Here,  AT  is  a temperature  change  in 

the  absence  of  plastic  relaxation. 

We  have  examined  three  different  cases  in  which  electron  microscopy 

M 

study  was  made  and  the  physical  properties  to  calculate  / are  available. 

These  are  (1)  A^O^  and  ZrO^  inclusions  in  Fe  matrix  [20],  (2)  a borosilicate 

glass  inclusion  in  AgCf  [21]  and  (3)  a carbide  inclusion  in  Mo  [22].  In 

the  last  example,  the  inclusion  was  not  positively  identified,  but  was  assumed 

M 

to  be  Mo0C.  Table  I lists  the  results  of  our  calculation  of  t together 

c ma  X 

with  experimental  data  and  comments  regarding  the  observations  of  dislocations. 

If  dislocation  loops  are  to  be  generated  at  the  critical  shear  strength 

in  a perfect  crystal  of  about  p/30  [23],  the  present  results  indicate  that  a 

borosilicate  glass  inclusion  in  AgC£  produces  dislocation  loops  around  it  and 

that  no  dislocation  is  produced  around  an  inclusion  in  Fe.  These  agree 

with  the  experiments.  However,  some  irregular  dislocations  were  observed 

M 

around  a ZrCL  inclusion  in  Fe.  The  calculated  result  show  that  t is  much 

£ Ilia  X 

less  than  p/30,  suggesting  that  irregular  dislocations  observed  around  ZrO^ 
inclusions  are  produced  by  non-thermal  stress.  As  for  the  third  example,  the 
calculation  indicates  that  M02C  inclusion  is  not  likely  to  create  dislocations. 
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This  suggests  either  that  prismatic  loops  observed  in  the  experiment  has 
resulted  from  the  interaction  with  glide  dislocations,  or  that  the  inclusion 


is  not  Mo^C. 


5.  INHOMOGENEITY  EFFECT 

If  a perfectly  fitting  inclusion  with  no  misfit  strain  has  elastic 
constants  differieng  form  those  of  the  matrix,  the  elastic  field  caused  by 
an  external  source  will  be  distrubed  around  the  inclusion.  Suppose  an  uni- 
axial external  stress  is  applied  along  X axis  in  Fig.  1.  The  applied  stress 
A A A 

o is  equal  to  and  other  terms  of  a.. . is  zero.  It  produces  elastic 

A 

strain,  e..,  which  is  uniform  at  large  distances  from  the  inclusion,  as 

* J 


(18) 


A A 

where  e is  given  by  a /E.  We  want  to  obtain  how  this  elastic  field  is 
perturbed  by  an  oblate  inclusion  which  has  different  elastic  constants  from 
those  of  matrix. 

5-1.  Internal  Stress  within  the  Inclusion,  (oi.).  L 
1 — 1 — i j -*-1  n h 

Eshelby  [8]  sol ved  the  inhomogeneity  effect  for  an  ellipsoidal  inclusion 
in  the  matrix  under  an  applied  stress.  He  showed  that  internal  stress  within 
the  inclusion,  (a|.)^n^,  can  be  written  as  Eq.  (19)  instead  of  Eq.  (2). 


^°i j^inh  " CijkiL(ekt  + ekt  ' ekJ^  ~ Cijk^ekt+  eki^- 


'19) 


From  Eqs.  (18)  and  (19),  ej.  for  the  inhomogeneity  effect  are  obtained  taking 

’ J 


12 


r 


v = v*  = 1/3  as. 


and 


eT  = eT  = x . = (1  - m)  {(1  - m)  r(k)  - 1}  . eA 

11  22  mls  3{(i  - m)2  fU)  : (1  - .)  fl(k)  ♦ » 

eT  =z.  (1  - m)  {3  - (1  - m)  s(k)}  CA 

33  inh  3{(1  - m)2  f(k)  - (1  - m)  g(k)  + 1} 

= o ifii'o 


(20) 


where  r(k)  - 3S1133  + S3333  and  s(k)  - 35^^  + 3S]122  + 2S33]1 

values  of  r and  s are  presented  in  Fig.  7. 

The  resultant  (aJ.).  . can  be  written  as 
ij  inh 

(°11  ^inh  = ^°22^inh  = Zu*  * (1  - m)  * ^inh  + zinh^ 

’ (l  - m)  * ^xinh  + Zinh^ 

0 if  i f j . 


and 


^°33^inh 


Numerical 


(21) 


By  using  x.  ^ and  z^nh  in  Eq.  (20),  the  internal  stresses  were  calculated 

as  a function  of  k.  The  results  are  plotted  against  k in  Figs.  8 and  9 for 

five  different  values  of  m.  The  stresses  were  normalized  by  o . When  m = 1, 

I I A 

one  obtains  (o^).^  = 0 and  (o33)inh  = o , which  coincide  with  the  applied 
stress.  (033)^  becomes  equal  to  oA  when  k vanishes  regardless  of  m.  As 
k approaches  unity,  the  effect  of  m on  (033)^  increases.  By  ploting  (a33)inh 
for  k > 1,  it  can  be  shown  that  the  stress  approaches  an  asymptotic  value 
depending  on  m at  k = «>.  On  the  other  hand,  (o|-|)inh  approaches  an  asymptotic 
value  as  k goes  to  zero  (see  Fig.  8).  Note  that  the  ratio  becomes 
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negative  when  the  inclusion  is  stiffer  than  the  matrix  (m  > 1).  This  ratio 
remains  positive  for  m < 1. 


I 


5-2.  Stress  Concentration  at  Matrix-Inclusion  Boundary 

When  we  consider  the  internal  stress  of  the  matrix  due  to  the  inhomo- 
geneity effect,  the  most  interesting  stress  compoent  at  the  boundary  is 
^°33'inh  at  ecluator  P°int  A-  This  becomes  equivalent  to  the  stress  intensity 
factor  problem  of  Mode  I crack  opening,  by  putting  m = 0.  At  equator  point  A, 
cose  in  Eq.  (5)  is  equal  to  1 , so  a^in  Eq.  (6)  becomes  as 


°33  + e33^  - 3x  - z}  . 

u 

Hence,  (033)^  at  A can  be  written  as, 

(°33^inh  = rE^{4(Sllll  + S1122  + 2S331 1 ^ " 3}xinh 

+ {4(S3333  + S1133)  " 1}zinh  + f‘e^  * ^22^ 


By  substituting  x..^  and  z^nh  from  Eq.  (19),  (0^3)^^  was  calculated  as  a 
function  of  k.  The  results  are  prented  in  Fig.  10  for  four  different  values 
of  m.  At  m = 0,  Eq.  (22.;  becomes  simply. 


, M , .A 
°33  ini/° 


r(k)  - 3s(k)  + 8 
8{1+  f(k)  - g( k) } 


(23) 


When  k is  zero,  we  have  f(0)  = 0,  and  g(0)  = r (0)  = s ( 0)  = 1 , so 

becomes  infinte.  At  k = 1 , we  have  (ol).  u/oA  = 2.0625  with  v = 1/3.  When 

oo  inn 

M A 

k becomes  smaller  than  0.01  at  m = 0,  can  be  approximated  by  1.28/k. 


i 


- 14  - 


The  radius  of  carvature  for  an  ellipse  in  Fig.  1 is  given  by 


{a4  + (c2  - a2)X2}' 


2 M 

At  equator  point  A(X1  = a),  p = c /a,  therefore  (033)^  for  a flat  cavity 
can  be  rewritten  in  terms  of  p as  = 1.28  x x iafp  . The  stress 

intensity  factor  of  "Mode  I"  crack  opening  is  given  as  [3], 


Kt  = Him  — 

1 p - 0 2 


By  substituting  (033)^*  Ki  f°r  a flat  cavity  can  be  written  as, 


1 , _ 1 . 28/n  A j — ~ t 1 t A iz 

Kj  = 2 'a  **  ~ 1.13  a /a 


This  agrees  with  Sneddon's  result  for  a penny-shaped  crack  [3,27].  The 
stress  distribution  around  a spheroidal  cavity  was  obtained  by  Neuber  [5] 
for  various  types  of  stressing.  In  the  case  of  simple  tension  in  the 

M A 

direction  of  the  axis  of  revolution  (a33)-jnt/a  vs-  a/P  were  plotted  using 
three  specific  values  of  v.  Our  corresponding  result  at  m = 0 (Fig.  8)  agrees 
well  with  his  results  taking  v = 0.3. 

The  evaluation  of  (033)^  for  an  inclusion  (m  f 0)  is  important  in 
relation  to  fatigue  and  brittle  fracture.  As  can  be  seen  in  Fig.  10, 

M 

(033)^  decreases  with  increases  in  the  value  of  m.  Its  dependence  on  m at 
a given  k becomes  quite  large  as  k decreases.  When  k is  equal  to  zero. 


can  be  written  as, 

33  inh 


°33  inh 
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M A 

Using  Eq.  (24),  values  of  m at  selected  ratios  of  were  obtained 

and  listed  below: 


m 

0 

0.07 

0.46 

1 

3 

10.3 

38 

00 

, M v , A 

(o33}inh/a 

CO 

10 

2 

1 

0 

—CO 

It  can  be  seen  that  the  stress  concentration  factor  becomes  less  than  2 
when  the  shear  modulus  ratio  m is  between  0.46  and  10.3.  Note  that  when  m 
is  greater  than  3,  the  sign  of  the  stress  ratio  becomes  negative;  i.e.,  when 
a compressional  stress  is  applied,  a large  tensile  stress  results  at  equator 
point  A for  a thin  disc  inclusion  having  a large  m value. 

Edwards  [6]  evaluated  the  stress  concentration  factor  of  prolate 
spheroidal  inclusion  (k  > 1).  Our  results  in  Fig.  10  agree  well  with  his 
at  k = 1 . 

6.  PLASTIC  DEFORMATION  EFFECT 

Suppose  a plastically  non-deformable  inclusion  exists  in  a deformable 
matrix.  When  applied  stress  deforms  the  matrix  plastically  internal  stresses 
develop  around  the  inclusion.  Tanaka  and  Mori  [9]  evaluated  such  internal 
stresses  and  developed  a work-hardening  theory  of  composite  materials  by 
using  Eshelby's  method.  However,  only  the  limiting  cases  were  considered. 

In  order  to  compliment  the  present  calculations,  it  is  necessary  to  evaluate 
effects  of  k for  quantitative  analysis  of  the  plastic  deformation  effect. 
Following  Tanaka  and  Mori  [9],  the  uniform  plastic  deformation  of  the  matrix 
proceeds  by  multiple  glide,  producing  uniaxial  strain  e along  x^  axis.  The 
geometry  of  an  inclusion  is  that  given  in  Fig.  1.  We  can  write  a stress-free 
transformation  strain  inside  the  inclusion  as  follows: 
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1/2  0 


0 -1 


Ashby  [28]  first  described  the  equivalent  condition  by  using  a 
dislocation  concept;  i.e.,  the  generation  of  dislocation  loops  at  the  inter- 
face of  the  inclusion.  In  the  method  of  Tanaka  and  Mori,  the  elastic  fields 
are  obtained  by  employing  Eq.  (2).  Note  that  as  a consequence  of  the  stress 
free  strain  [Eq.  (25)],  a spheroidal  inclusion  with  the  aspect  ratio  k is  to 
be  considered  as  the  one  with  k'  = k(l  + 2 • ep)  Tor  the  first  approximation  [9], 


6-1.  Internal  Stress  Within  the  Inclusi 


125^4  A 


e!  . in  Eq.  (2)  for  plastic  deformation  effect  can  be  expressed  as 

* J 


el  = eT 
*11  2 


,T  _ x _ m{l  - (1  - m)u(k) } . £ 

22  " p ’ 2{(1  - m)2f(k)  - (1  - m)  g(k)  + 1}  p 

m{(l  - m)  v ( k)  - 1} . r 


{(1  - mr  f(k)  - (1  - m)  g(k)  + 1} 


where  u(k)  - 2S1133  + S3333  and  v(k)  + S1122  + S3311  (see  Fig.  7). 

Internal  stress  within  the  inclusion,  (ol.)  , can  be  calculated.  The 

■ J P 

results  indicate  that  (a!.)  = 0 if  i / j,  and  the  diagonal  stress  components 

1 J r 


<4i>p  ' ■ 2v'  • 


•«3*p  * V ■ Ep/2) 


l^p  ■ • 


l(Xp  + *p>  * V2> 
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Numerical  results  are  presented  in  Fig.  IV  and  Fig.  12^,  where  the 

stresses  are  normalized  by  the  factor  E**ep.  As  k goes  to  0 (disc),  each 

of  (°ji)p/£**ep  ancl  (033)p/E**ep  approaches  to  an  asymptotic  value;  -0.75  and 

0,  respectively.  As  k goes  to  1 (sphere)  at  constant  m,  the  magnitude  of 

(on)n  decreases,  but  that  of  (oL)  increases.  At  k = 1 , we  have  (oL)  = 

1 1 P 33  p 33  p 

T 

2(o^)p  for  a given  m [10]. 

M 

6-2  Internal  Stress  at  the  Boundary,  (ov .) 

■ J P 

As  the  amount  of  plastic  deformation  of  the  matrix  is  increased,  the 
internal  stress  at  the  matrix-inclusion  boundary  develops  according  to  Eq.  (6). 

M 

For  example,  an  internal  stress  component,  (a33)p,  at  equator  point  A can  be 
written  as 


^33^  = 8 Et{4(Sllll  + S1 122  + ZS3311^  ' 3}xp 


+ {4^S3333  + S11 33 ^ ' 1}zp^ 


The  identical  component  at  polar  point  B is  given  by 


(a33}p  = I • E • • {<xp  + zp}  + V2}  • 

Figure  13  shows  the  results  obtained  by  substituting  xp  and  zp  from  Eq.  (26). 
At  polar  point  B,  we  have  o*3  = o^.  At  equator  point  A,  (0^3 ) p/E  • ep 
depends  on  the  values  of  m and  k.  When  m _>  1 , the  maximum  of  (o33)p  occurs 
at  k = 0.  In  this  case,  the  thinner  is  the  inclusion,  the  smaller  amount  of 
cp  will  initiate  the  decohesion  at  the  boundary.  On  the  other  hand,  when 

^The  horizontal  axis  In  the  figures  Is  to  be  k'  = k(l  + 3 cp/2),  however,  if 
e is  less  than  6X,  the  difference  between  k'  and  k becomes  less  than  10%. 

P 


* 


m < 1,  the  maximum  stress  is  reached  with  0 < k < 1;  for  example,  the 

maximum  of  (a,0)  occurs  at  k=0.2  when  m = 0.5. 

33  p 

As  in  the  cases  discussed  in  Sections  4 and  5,  not  all  the  non-diagonal 
/ M , 

components  of  (o. .)  vanish.  Thus,  shear  stresses  at  the  boundary  can  be 

* J r 

deduced,  although  extended  computations  are  required.  In  the  case  of  a 
spherical  inclusion,  Mori  and  coworkers  [29]  determined  the  shear  stresses 
and  their  maxima  at  the  boundary.  They  employed  the  calculation  in  the 
analysis  of  plastic  relaxation  around  a particle. 


7.  DISCUSSION 


The  present  study  has  evaluated  the  internal  stresses  within  inclusion, 

I M 

at.,  and  those  at  the  matrix-inclusion  boundary,  o': . , due  to  the  misfit,  in- 
1 J 1 J 

homogeneity,  and  plastic  deformation  effects.  Specifically,  = o*,,,  033 

M 

and  a33  have  been  obtained  as  a function  of  k and  m for  an  oblate  spheroidal 

inclusion.  Now,  one  can  obtain  the  internal  stresses  once  the  values  of  m and 

k are  known.  While  the  orientation  of  the  inclusion  has  no  essential  effect 

on  the  misfit  induced  stresses  because  of  isotropy,  the  inhomogeneity  and 

plastic  deformation  effects  are  limited  to  the  specific  geometry  considered. 

/\ 

In  the  present  study,  the  direction  of  external  stress  o coincides  the  axis 
of  revolution  of  the  inclusion,  which  is  also  the  direction  of  uniaxial 
plastic  strain  ep. 

Within  the  limits  of  the  present  calculation,  a variety  of  problems  in 
materials  research  can  be  analyzed  in  addition  to  those  examined  in  the  previous 
sections  of  this  paper.  Since  mechanical  responses  of  inclusions  in  structural 
materials  are  of  significant  concern,  we  consider  as  an  example  the  internal 
stresses  in  and  around  an  MnS  inclusion  in  steel. 
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It  is  well  known  that  the  decohesion  of  flattened  oblong-shaped  MnS 

inlcusions  in  steel  plates  is  responsible  for  the  loss  of  ductility  along 

the  thickness  direction.  The  loss  is  particularly  pronounced  when  the 

recently  developed  controlled  rolling  process  is  adopted.  Acoustic  emission 

studies  of  high  strength  low  alloy  steels  have  revealed  that  the  decohesion 

or  fracture  of  the  inclusions  takes  place  well  before  the  yield  stress  is 

reached.  [30,31].  Much  of  decohesion  are  complete  within  the  initial 

several  percent  of  plastic  deformation.  Let  us  consider  the  stress  in  and 

around  an  MnS  inclusion.  Its  shape  is  approximated  by  an  oblate  spheroid  and 

the  applied  stress  is  along  the  axis  of  revolution.  It  is  assumed  that  the 

matrix  ceases  to  allow  plastic  relaxation  around  the  inclusion  below  800°K, 

giving  rise  to  thermal  misfit  strain.  The  thermal  expansion  coefficients  of 

MnS  and  steel  are  found  to  be  18.1  x 10~*7°K  and  12.0  x 10~*V°K,  respectively 
T _ q 

[32],  which  gives  e = -3  x 10  . Young's  moduli  of  MnS  and  steel  are  137 

GPa  and  206  GPa,  respectively,  which  correspond  to  m = 2/3  [32].  The  magnitude  of 
was  taken  as  E/1000  or  206  MPa.  The  amount  of  Cp  was  set  at  0.1%  or  1%, 

A 

although  o would  be  higher  than  E/1000  when  the  latter  value  was  employed. 

The  internal  stress  o*3  inside  the  MnS  inclusion  was  calculated  for  the 
three  different  effects  (cf.  Eqs.  8,  21  and  27).  The  results  are  presented 
against  k in  Fig.  14.  Individual  components  as  well  as  the  sum  of  the  three, 

<°33>tot*l  = (c,33)mis  + (033>1nh  + <°33V  are  9iven'  As  k at,proaches  2er°- 
the  contributions  of  (a33)miS  and  ^33^  become  negligible.  Plastic 

deformation  effect  is  small  regardless  of  k until  cp  exceeds  several  percent. 

Misfit  effect  becomes  significant  as  k increases;  (°33)m-js  exceeds  (°33)-jn^  for 

K > 0.18.  For  a typical  aspect  ratio  of  MnS  inclusions  of  about  0.1,  both 

misfit  and  inhomogeneity  effects  produce  the  normal  stress  on  the  broad  face 

of  the  inclusion.  The  internal  stress  0^3  at  equator  point  A of  the  matrix- 

inclusion  boundary  was  also  obtained  by  using  Eqs.  (6)  and  (22).  The  results 
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are  shown  in  Fig.  15.  As  can  be  seen  in  Fig.  15,  the  total  internal  stress 

at  point  A,  (oS3)total  = (o^mis  + (o33}inh  + (o33}p’  is  ne93tive  as  far  as 

e is  less  than  2%.  This  is  due  to  a large  negative  values  of  the  misfit 
M 

effect.  If  (a,J  . remains  unrelaxed  by  plastic  deformation  of  the  matrix, 
it  is  most  unlikely  that  the  decohesion  or  fracture  begins  at  the  edge  of  an 
inclusion.  Even  if  some  plastic  relaxation  occurs,  the  present  analysis 
indicates  that  the  normal  stress  on  the  inclusion  plane,  (033)^  and  (033)^5 
in  particular,  contributes  most  significantly  to  the  decohesion  of  the  inclusion. 

Tanaka  et  al . [10,11]  were  the  first  to  apply  Eshelby  theory  to  the 
decohesion  problem.  However,  misfit  effect  was  not  included.  When  a 
spherical  inclusion  is  considered,  their  conclusion  is  basically  valid. 

Because  of  the  omission  of  the  misfit  effect,  the  contribution  of  plastic 
deformation  effect  appears  to  be  overemphasized.  In  dealing  with  the  disc- 
shaped inclusions,  however,  the  geometry  was  inappropriate  in  that  the  deco- 
hesion was  considered  to  occur  at  the  edge  of  the  an  infinitely  thin  inclusion 
via  oj-..  Again,  the  plastic  deformation  effect  was  overestimated.  Thus,  their 
conclusion  concerning  disc-shaped  inclusions  is  inapplicable  to  configurations 
of  general  interest. 
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FIGURE  CAPTIONS 


Fig.  1.  The  orientation  of  an  oblate  spheroidal  inclusion.  Applied  stress, 

A 2 2 2 

a , is  exerted  along  the  X^.  Equator  point  A refers  to  + X-j  = a 

and  polar  point  B refers  to  (0,0,c). 

Fig.  2.  Shape  dependence  of  parametric  functions;  f,  g,  h and 
Fig.  3.  Normal  components  of  internal  stress  (normalized  by  E c ) inside  the 
inclusion  due  to  misfit  effect  vs.  aspect  ratio,  k.  Tne  ratio  of 
Young's  moduli  employed  is  indicated. 

Fig.  4.  Maximum  shear  stress  inside  the  inclusion  due  to  misfit  effect. 

^iS/u*eT  against  k. 

Fig.  5.  Elastic  strain  energy  of  the  inclusion  and  matrix  due  to  misfit  effect. 

E /4p  (£T)2Vt  against  k. 

S I 

Fig.  6.  Internal  stress  at  the  matrix-inclusion  boundary  due  to  misfit  effect. 

a9ainst  k- 

Fig.  7.  Shape  dependence  of  parametric  functions;  r,  s,  u and  v. 

Fig.  8.  Internal  stress  inside  the  inclusion  due  to  inhomogeneity  effect. 

(°11 )inh/°A  dgainst  k- 

Fig.  9.  Internal  stress  inside  the  inclusion  due  to  inhomogeneity  effect. 
(033)inh/°A  a9ainst  k- 

Fig.  10.  Internal  stress  at  the  equator  point  A on  the  matrix-inclusion 

M A 

boundary  due  to  inhomogeneity  effect.  ^°33  1nl/°  a9a,'"st  k- 

Fig.  11.  Internal  stress  inside  the  inclusion  due  to  uniaxial  plastic  deformation 

I * 

of  the  matrix.  (o^)p/E  ep  against  k. 

Fig.  12.  Internal  stress  inside  the  inclusion  due  to  plastic  deformation  effect. 

( °33 ) p/E  ep  a9ainst  k. 

Fig.  13.  Internal  stress  at  equator  point  A and  polar  point  B at  the  matrix- 

inclusion  boundary  due  to  plastic  deformation  effect.  °33/^f-p  against  k. 


Fig.  14.  Internal  stress  inside  an  oblate  MnS  inclusion  in  the  steel  matrix 

(m  = 2/3).  Three  components  of  and  their  sum  for  the  case  of 

A -3  -3  * 

Ep  = 1%  are  shown  against  k.  a was  taken  as  10  E or  1.5  x 10  E 

Fig.  15.  Internal  stress  at  equator  point  A of  the  matrix-inclusion  boundary 
M 

o^/E  against  k.  Note  a large  negative  contribution  due  to  misfit 


effect. 


Fig.  1.  The  orientation  of  an  oblate  spheroidal  inclusion.  Applied  stress, 

2 2 

oA,  is  exerted  along  the  X3-  Equator  point  A refers  to  X1  + = a 

and  polar  point  B refers  to  (0,0,c). 


Fig.  3.  Normal  components  of  internal  stress  (normalized  by  E cT)  inside  the 
inclusion  due  to  misfit  effect  vs.  aspect  ratio,  k.  The  ratio  of 
Young's  moduli  employed  is  indicated. 
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Fig.  4.  Maximum  shear  stress  inside  the  inclusion  due  to  misfit  effect. 
I * T 

Tmis^  e a9ainst  k. 


Elastic  strain  energy  of  the  inclusion  and  matrix  due  to  misfit  effect 
E_/4 u*(eT)ZVT  against  k . 
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Fig.  6.  Internal  stress  at  the  matrix-inclusion  boundary  due  to  misfit  effect. 
(o33)mis/EeT  against  k- 
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Internal  stress  inside  the  inclusion  due  to  inhomogeneity  effect 
(°!i)<nh/oA  against  k . 


Fig.  10.  Internal  stress  at  the  equator  point  A on  the  matrix-inclusion 

M A 

boundary  due  to  inhomogeneity  effect,  against 
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Fig.  11.  Internal  stress  inside  the  inclusion  due  to  uniaxial  plastic  deformation 
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of  the  matrix.  (oj-|)p/E  ep  against  k. 
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Fig.  12.  Internal  stress  inside  the  inclusion  due  to  plastic  deformation 
effect.  (°33)p/E*  ep  against  k. 
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Fig.  13.  Internal  stress  at  equator  point  A and  polar  point  B at  the  matrix- 

inclusion  boundary  due  to  plastic  deformation  effect.  °33^cp  gainst  k. 
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Fig.  14.  Internal  stress  inside  an  oblate  MnS  inclusion  in  the  steel  matrix 

(m  = 2/3).  Three  components  of  and  their  sum  for  the  case  of 

A .3  ★ 

Cp  = 1%  are  shown  against  k.  0 was  taken  as  10  E or  1.5  * 10  E . 


L 


